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VIEWING ANGLE AND ENVIRONMENT EFFECTS IN GRB: 
SOURCES OF AFTERGLOW DIVERSITY 

P. Meszaros \ M.J. Rees^ & R.A.M.J. Wijers^ 
§^ ■ ABSTRACT 

^ I We discuss the afterglows from the evolution of both spherical and anisotropic 

Q ' fireballs decelerating in an inhomogeneous external medium. We consider both the 

• radiative and adiabatic evolution regimes, and analyze the physical conditions under 

which these regimes can be used. Afterglows may be expected to differ widely among 
^ . themselves, depending on the angular anisotropy of the fireball and the properties 

, of the environment. They may be entirely absent, or may be detected without a 

, corresponding 7-ray event. A tabulation of different representative light curves is 

0^ I presented, covering a wide range of behaviors that resemble what is currently observed 

^ ■ in GRB 970228, GRB 970508 and other objects. 



Subject headings: gamma-rays: bursts 

o 
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^ ■ 1. Introduction 

^ I The discovery of the afterglows of gamma-ray bursts (GRB) provides new information which 

H ' can constrain the models used to explain these objects. Significant interest was aroused by the fact 

that several of the features reported in the first GRB detected over time scales <> days at X-ray (X) 



and optical (O) wavelengths, GRB 970228 ( posta et al, 199"7a ) agreed quite well with theoretical 



expectations from the simplest relativistic fireball afterglow models published in advance of the 
observations ( [Meszaros fc Rees, 1997a| ; see also Vietri, 1997a| ). A number of theoretical papers 



were stimulated by this and subsequent observations (e.g. Tavani, 1997; Waxman, 1997a;Reichart, 



1997; Wijers, Rees & Meszaros , 1997, among others), and interest continued to grow as new 



observations provided controversial evidence for the distance scale and the possible host ( |Sahu et 



al, 1997| ). New evidence and new puzzles were added when the optical counterpart to the second 



discovered afterglow (GRB 970508) was attributed a cosmological redshift ([Metzger et al., 1997|) , 
as well as a radio counterpart ( [Frail, et al, 1997 ; p?aylor, et al, 1997 ) and X/0 light curves showing 



^Dpt. of Astronomy & Astrophysics, Pennsylvania State University, University Park, PA 16803 
^Institute of Astronomy, Cambridge University, Madingley Road, Cambridge CBS OHA, U.K. 



- 2 - 



a rise and decline ( Djorgovski, S. et al., 1997| ; |Fruchter, et al, 1997] ). Some bursts, however, were 
detected only in X but not O (e.g. GRB 970828), while some which would have been expected to 
be seen in X or O were not (e.g. GRB 970111). 

Additional structure on the light curves has also emerged from a continued analysis of some 
of these objects down to the faintest flux levels. The large variety of behaviors exhibited by 
afterglows, while clearly compatible with relativistic fireball models, poses new challenges of 
interpretation, e.g. Waxman, 1997b| ; Vietri, 1997b ; Katz &: Piran, 1997 ; Rhoads, 1997 ; Paczyhski, 



1997 



Some of the questions at the forefront of attention include the effect of the external 
medium, the degree to which afterglows may be considered to be isotropic events, and the effects 
of the radiative efficiency on the evolution of the remnant. We address all three of these issues 
here. We also clarify some of the issues that have been recently raised about the dynamical 
effects of different radiative efficiency regimes. We then discuss the possible variety of afterglow 
behavior that is expected from isotropic or anisotropic fireballs expanding in a medium which 
may be inhomogeneous, either due to external gradients or due to expansion in an irregular 
cavity. We apply these models to interpret some of the salient observational features of several 
GRB afterglows, and discuss their possible use for predicting detection rates of X/O/R afterglows 
undetected in 7-rays, as well as some possible reasons for the non-detection of afterglows in GRB. 



2. Expansion Dynamics and Radiative Efficiency 

In some bursts (e.g. GRB 970508) the afterglow seems to contain a significant amount of 



energy compared to the typical (isotropic) estimate of ~ 10^^ erg s ^. This led Vietri, 1997b 



to suggest that the afterglow must remain radiatively efficient ~ weeks after the burst and evolve 



with r oc r~'^. This regime was also considered in [Katz Sz Piran, 1997| , who refer to previous 
relativistic fireball models as radiating only a small fraction of the total kinetic energy, and go on 
to consider instantaneously cooling fireballs. It is important to discuss in more detail what are the 
conditions necessary for the radiative efficiency having an effect on the dynamic evolution of an 
expanding cloud. 

The classical fireball models, with "isotropic equivalent" energies of ~ 10^^ erg and bulk 
Lorentz factors F ~ 10^ — 10^ have, in fact, been generally been taken to be in the radiative stage 
during the 7-ray event, i.e. radiative efficiency near unity, meaning that of order the initial total 



kinetic energy of the protons is radiated in the observer- frame expansion time (Rees & Meszaros 



1992 , and subsequent papers). However, for some parameters the bulk of this energy can appear 



at energies other than MeV ([Meszaros , Rees fc Papathanassiou, 1994 ; Meszaros &: Rees, 1994 ). 



This high efficiency in the initial deceleration shock can occur if the electrons are assumed in the 
shock to be heated to 7e ~ ^e{'nip/me)V with reasonably close to unity. Experimental evidence 
from interplanetary collisionless shocks indicates that this could be the case. In such fireballs 
the electrons are likely to retain high radiative efficiency for some time after the GRB, and 
energetically the most important are the newly shocked electrons near the downwards evolving 
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peak, of initial post-shock energy 7e ~ ^e("i-p/'TT'e)r(t). This peak is generally where most of 
the electron internal energy is. In the regime where the peak electrons have high radiative 
efficiency, if throughout the entire remnant volume the protons are able to establish ( and remain 
in) equipartition with the electrons, then the remnant evolves with F oc r~^ in a homogeneous 



external medium (e.g. Blandford & McKee, 1976; Vietri, 1997a; Katz & Piran, 1997). This follows 



simply from the momentum conservation law, if we can assume that the radiative losses tap also 
the proton and magnetic energy (strong coupling), and the radiative time scale is shorter than the 
dynamic time scale, so that energy conservation cannot be used. This is the classical "snowplow" 
approximation of supernova remnants. The alternative regime is that where the radiative losses 
do not tap the proton and magnetic energy, only the electron energy (weak coupling), and/or the 
radiative cooling time scale is longer than the dynamic time scale. In this case, possibly after an 
initial short cooling of the electrons, one can assume energy conservation (most of the energy is 
in the protons and/or magnetic fields), and one has T oc r~^/^ in a homogeneous medium (e.g. 
Blandford fc McKee, 1976t [Paczyiiski fc Rhoads, 1993| , |Katz, 1994b| , pVIeszaros fc Rees, 1997aD . 



The (strong coupling) radiative regime and the adiabatic regime can be generalized to the 
case where the fireball moves into an inhomogeneous external medium. We consider a spherical 
fireball of energy E and bulk Lorentz factor T which are independent of angle 6, expanding into an 
external medium of density n(r) oc r~'^ which is also independent of the angle, where r is distance 
from the center of the explosion. The shocked gas evolves according to the conservation law 

pi+A^3^ OC T^+\^-<^ oc constant , (1) 

where A = 1(0) corresponds to energy (momentum) conservation, i.e. to the adiabatic (radiative) 
regimes. These regimes must be understood in a global or dynamic sense, as applying to the entire 
remnant, i.e. baryons, magnetic fields, electrons, etc., or at any rate to its dynamically dominant 
constituents. Since the observer- frame or detector time t must satisfy r oc ciT^, we have 

r oc r--(3-'i)/(l+^) Of ^-{■i-d)/(J+A-2d) ^ ^ ^ ^(l+A)/(7+A-2d) 1^2) 

If (i = and a remnant starts out in the strong coupling radiative regime {A = 0) then F oc r~^. 
However after the expansion has proceeded for some time eventually the cooling time of the 
electrons at the peak of the distribution becomes longer than the expansion time, at which 
point energy conservation (adiabatic approximation) becomes valid (eq.|l| with ^ = 1) leading to 
F oc r~^/^. The power law of equation (||) is valid as long as the remnant is relativistic, i.e., until 
the time tnr/to ^ '^^o'^^ 2d)/(3 d) ^ where to ~ t^ is the duration of the 7-ray burst itself and Fq is 

.{7+A)/3 



the initial Lorentz during the burst. For a homogeneous medium d = and tnr/^o ~ Fy^"^^^ e.g. 



for an adiabatic remnant A = \ with Fq = 10^ and ^ = 1 s, one has tnr ~ 1 year. The subsequent 



behavior in the nonrelativistic stage is described in Wijers, Rees &; Meszaros , 1997. Shorter times 
for reaching the nonrelativistic stage are possible for a radiative remnant, or for a radiative stage 
followed by an adiabatic one, while longer times can occur for longer or for expansion into a 
medium whose density decreases with distance (0 < d < 3). 
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However, the conditions under which the remnant dynamics may be considered adiabatic or 
radiative is far from unambiguous, and is crucially dependent on poorly known questions about 
post-shock energy exchange between protons and electrons. Electrons cool quickly compared to 
protons, and it remains an unsolved question, of importance also in other areas of astrophysics, 
whether behind the shocks, after the electrons have cooled, the protons remain hot (i.e. a two 
temperature plasma, as in hot torus models of AGN), or whether they tend toward some degree 
of equipartition with the cooled electrons by virtue of unknown fast energy exchange mechanisms. 
The collisionless shock transition is the only place where one is guaranteed fast varying chaotic 
electric and magnetic fields which can lead to quick relaxation, and it is known that an interplay 
between protons and magnetic fields can occur there which can lead to values close to equipartition. 
Such equipartition is also inferred from measurements in the ISM. However, in fast evolving flows 
such as those in GRB, it is unclear whether such equipartition can occur anywhere except possibly 
near the shock transition itself. If behind the shocks the protons are unable to quickly readjust to 
the electron losses, then even if the electrons are radiative (a = 1) the protons can be adiabatic 
(A = 1) leading to F oc r~^/^, in a homogeneous external medium. A related problem is that, in 
order for the remnant to evolve with T cx (again for a homogeneous external medium), the 
magnetic energy should also decrease, since if the latter were conserved it would soon dominate 
the total energy density and the remnant would evolve as a polytrope with adiabatic index 4/3 
which leads to F oc r"'^/^. Thus both the proton and the magnetic energy need to be transferred 
on a fast time scale to the electrons in order to ensure F oc r~^. An additional factor that might 
contribute towards a steepening of the decay of F are other energy losses, e.g. such as from the 
escape of accelerated nonthermal protons from the shell, if these carry substantial energy. So far, 
there are neither detailed simulations nor experimental evidence concerning this in GRB. In the 
absence of such losses, or of a quick energy exchange between post-shock protons plus magnetic 
fields and the electrons, one can therefore have a situation where the electrons are "radiatively 
efficient", but the dynamics of the remnant expansion follows an "adiabatic" law. This occurs 
if the electron cooling time is less than the expansion time. The shocked electrons can radiate 
up to half of the proton energy in the shocks, but the protons and the magnetic fields retain at 
least half and this would be enough to ensure a quasi-adiabatic dynamic evolution of the remnant 
with F oc r~'^/^. The latter is also true when the electrons are adiabatic. For an inhomogeneous 
medium that decays with radius, both of these decay laws would be flatter. 

There is no difficulty in treating the weak coupling case where the electrons are radiatively 
efficient but the dynamics is adiabatic (§^). This regime is physically as plausible, if not more, as 
the strong coupling one where protons and fields exchange energy with electrons on a fast time 
scale. For reasons of simplicity, in |^ and in the rest of the paper we will assume that magnetic 
fields are near equipartition with the protons, which ensures a simple expression for the electron 
radiative efficiency in terms of only the synchrotron cooling time and the expansion time (the 
situation where inverse Compton (IC) cooling is important only introduces some extra changes in 
the way the synchrotron efficiency is defined). 
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3. Spherical Inhomogeneous Models 

As in the previous section, we consider a spherical fireball of energy E and bulk Lorentz 
factor r independent of angle d expanding into an external medium of density n(r) oc r^'^. In 
the simplest afterglow model one considers the time evolution of the radiation from the external 
medium shocked by the blast wave as it slows down. Denoting quantities in the comoving frame 
of the shocked fluid with primes, in the post-shock region the density is n' ~ nV, the mean proton 
and electron random Lorentz factors are 7p ~ F and 7e ~ ^eC'T^p/'^e)!^ (where (me/mp) ^ ^ 1 
is the fraction of the electron equipartition energy relative to protons), the turbulently generated 
magnetic field (assumed to build up to a fraction .^_b < 1 of the field in equipartition with the 
random proton energy) is B' oc .^sn^/^r, and the peak of the electron synchrotron spectrum is 
at comoving frequency u'^ oc B'^y"^ oc ^B^gn^'^/^r^ oc r~'^/^r^. The corresponding observer peak 
frequency is 

OC e^dn'/'r^ oc i-[12-(^/2){7-A)]/(7+A-2d) (3) 

The synchrotron radiative efficiency at is esy,m ~ t'~y „^/ {t'^y „^ + t'^^ + t'^^^^j.) where 
^'sy,m ^ l/lCeC^r^n) is the comoving synchrotron cooling time at i/^, t'^^ is comoving expansion 
time or adiabatic cooling time and t'o^/jgr other loss mechanism, e.g. inverse Compton (IC), 

if important. In the limit where only synchrotron and/or adiabatic losses are important we may 
write esy^m ~ {^'ex/'^'sy,mY i which is unity in the electron radiative (a = 0) regime, and < 1 in the 
electron adiabatic (a = 1) regime. We have egy^m ~ (^B^e^'^-r^)'^ oc {i%ie'r^~'^^'^)°' ■ 

We first assume (§|2|) that protons and magnetic fields are strongly coupled to 
electrons, so if the electrons are radiative the entire remnant is radiative, and the 
index A = a. The comoving synchrotron intensity at the comoving peak frequency 

is II,^ oc KiPly/u'Jct'^,^ ~ n',{^,m,C^/u'Jcesy,m oc ^'^--^^a-l^a-d[a+l/2)^2a-l ^ ^^^^^ 

^'min ~ i'sy^sy,m is the shortest of the possible cooling times (synchrotron or adiabatic, in the above 
approximation). The flux from the relativistically expanding source at observer frequency is 
F,„cxt2r5/^^,or 

JP^^ oc z,^-[2(l-")-W2)(l+a)]/[12-{rf/2){7-a)] ^ ^[2{l-a)-{d/2){l+a)]/{7+a-2<i) ^ ^^-^ 

scaling with Cr^Cs"^- If the expansion is in the radiative a = regime, 
^i^ni 0^ ^rn^ ('^/2)]/[i2 (7/2)d] ^ ^{A-d) / (lA-Ad) increases in time for any d < 3, being oc Um^^^ oc t'^^''' 
in a homogeneous medium with d = (Vietri 1997b obtains a different scaling by taking in /^^ 
the shocked gas comoving width AR' as path length, which however is equal to ct' ■ only for 



adiabatic expansion, e.g. Meszaros &: Rees, 1997a , where t'^in ~ t'ex)- the adiabatic a = 1 



regime F^^ oc Vm^^'^ ^"^^ oc t '^^^^ which is a constant independent of Vm and of time for 
adiabatic expansion in a homogeneous d = Q medium ( [Meszaros &: Rees, 1997a ; Katz, 1994b ) 



but Fy^ decreases in time for an inhomogeneous medium with < d < 3 (for d > 3 the fireball 
encounters most of the external mass near its initial radius). For a power-law spectrum F^, oc z^", 
the flux at a fixed detector frequency vd is 

FD = F,^{vD/ymr OC i[2(l-a)+12a-(d/2)({l+a}+a{7-a})]/[7+a-2d] _ (5) 
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scaling as ^-i+'^-2a^^(i+a)+2a^ j.q^ ^ ^g^^jj £qj, ^j^g strong electron-proton coupling regime. 
For a typical synchrotron spectrum a ~ l/3(— 1) below(above) the break frequency I'mit), as the 
latter decreases in time the flux from radiative a = models in the detector frequency band vjj 
at times for which < Um ,1^0 > is i^z? oc t^/"^ , t~^'^^'^ in a homogeneous d = medium, and 

in an inhomogeneous d = 2 medium. Adiabatic a = 1 models give 
in a homogeneous medium, and F^ oc , in an inhomogeneous d = 2 medium. Other values 
can be calculated from eq. @ for different a before and after the break. Tables 1 and 2 give for 
the isotropic strong coupling case several examples in the next to last column. 

A different regime is obtained if one assumes that the protons and magnetic fields are not 
strongly coupled to the electrons behind the shocks. In this case the dynamics of the remnant as 
a whole is controlled by the index A in eqs. (|^,|2|), and as long as F ^ 1 the remnant is adiabatic 
with A = 1 and T oc r~^/^ for a homogeneous medium, whether the electrons are radiative or not, 
i.e. independent of a (see More generally, in this case T oc i-{3-c()/(8-2d)^ ^ ^ ^2/(8-2d^^ 

OC iBil t-(12-3.)/(8-2.) (g) 

In the synchrotron efficiency one must keep a separate index a^Ato account for the electrons 
being radiative or not. We have then V oc ^^-^^'^-^r''^'^-d)-dl'iY'^o--^ ^ and F^^ ~ ^^^^K^ is 

In a given fixed detector band vd one observes for a typical spectral shape oc z/" a 
time-dependent flux Fd oc Fy^ , or 

Fd oc ^g'^-l-")^«-l-2" ^[4(l-a)-d+a(12-3d)]/{8-2d) ^^-^ 

This is valid in the relativistic expansion regime. After a remnant becomes nonrelativistic 
(see below equation 2), the flux in a homogeneous d = 0, a = 1 medium would steepen to 
Fd oc t{3+i5a)/5 ^ t-12/5 for a = -1 ( iWijers, Rees fc Meszaros , 1997| ). The variety of time 



behaviors possible for relativistic isotropic models in both the strong and weak coupling cases 
depending on whether they are radiative or adiabatic is shown in the last two columns of Tables 1 
and 2. 

The time behavior given by eqs. (^,§) can be complicated by at least two effects. One is 
that, unless observations start after the peak electrons are adiabatic, at some subsequent time 
the value of a in eqs. (Q,|5|,^,^) switches from to 1 as the peak electrons become adiabatic. The 
second is that if the detector frequency ud > z^mithe flow (being controlled by particles radiating 
at I'm) can already be adiabatic, while the smaller number of higher energy particles radiating at 
I'D may still radiate efficiently. The frequency and time power law dependences of the flux before 
and after electron radiative inefficiency occurs are different. For a steady injection of accelerated 
electrons, the self-consistent electron energy power law index p in the presence of fast synchrotron 
losses is one power steeper than the injected spectrum, and the self-consistent synchrotron power 
law index a = (p — l)/2 is a half power steeper than for the adiabatic (negligible loss) case. If 
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the lowest energy electrons near the peak 7e ~ ^e(mp/me)r are radiatively efficient, all electrons 
above that are as well. As the remnant evolves, the first electrons to become inefficient are the 
lowest energy ones (in the peak corresponding to i^m(^))> and a flattening break by 1/2 power in 
the photon spectrum at frequency ^^(i) > i'm(^) moves to frequencies increasingly higher than 
Vmit)- The electron Lorentz factor at which the synchrotron time just equals the expansion time 
r/cT is 7;, OC r^^^'^^^^'^^l/^^"'"^) and the corresponding "adiabatic" photon frequency is 

n OC B'jIT OC nV2r272 r'^^-im ^ ^-[2-(3/2)d](l+A)/(7+A-d) ^ 

which can either decrease in time for d < 4/3 (including a homogeneous medium with d = 0) or 
increase for d > 4/3 (although it always increases with respect to fm(i))- For an external medium 
whose density drops with radius faster than d ^ 4/3, if initially Vb > t^d it will always remain so, 
and the spectral index remains adiabatic without change (until a much higher cutoff is reached 
where the acceleration becomes inefficient and the spectrum drops off exponentially). However for 
a homogeneous medium or one with d < 4/3, if initially vi, > ud the photon spectral index will 
at some later time steepen by 1/2 as sweeps through the observing band vo and the observed 
spectrum transitions from the adiabatic to the radiative regime. 



4. Anisotropic Inhomogeneous Models 

The observed afterglow temporal decays are conventionally fitted by power-laws, and it is 
interesting to explore how the decay slopes would depend on the angular dependence of the 
dynamically relevant quantities of a fireball. To that effect, we consider anisotropic relativistic 
outflows where both the energy per unit solid angle and the bulk Lorentz factor depend on the 
angle 6 as power-laws (at least over some range of angles), and also consider the external density 
distribution to depend on radius as a power-law, 

EoqO-^ , T o^e~^ , n OC r"'^ . (10) 

If there is a well defined jet, the normalizations of E and T may be different for material inside 
and outside the jet opening angle Oq- At each angle the outflow starts converting a significant 
fraction of its bulk kinetic energy into radiation when an external blast wave develops at the 
angle-dependent deceleration radius oc [£'/n(rrf)r^]^/^, at an angle-dependent observer-frame 
(detector) time t ^ r /cV'^. The 0-dependence of eq.([lO|) implies that the deceleration blast wave 
at different angles occurs at 

p ^ ^-k{3-d)/{8k-j-2dk) ^ ^ ^ ^{2k-j)/{8k-j-2dk) 
E OC t-^i^-d)/{8k-j-2dk) ^ Q ^ ^(3-d)/i8k-j-2dk) _ 

Depending on the normalization of ([To| ) and causality considerations, the radiation from the blast 
waves occurring at increasing 9 at successive times t can dominate the afterglow evolution (as 
opposed to the decay of E and T along the same as a function of time). For instance, if the event 
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has been detected at 7-rays, and there is a jet of opening angle do, the observer is presumably 
within angles ^ 60 from the axis. For T = To{0/9o)^^ , in order for subsequent blast waves at 
r = rd from ^ > to be observed at times t > to one needs F"^ ^ ^ to be satisfied, that is 
0/eo> (0oFo)V(fc-i), or 

t/to ^ (0„ro)(S'=-^-2*)/I(3-'^)('="l)l . (12) 

For values of A; < 1 the blast waves are detectable at all angles, but for k > 1 there are ranges of 
k,j for which ^oFq is limited to values ^ 5 — 10 in order to detect the blast wave at reasonable 
t/to- In this case the initial part of the afterglow may be due to the evolution in time of the 
gas responsible for the burst initially observed, until such a time when the causality condition 
is satisfied for gas at larger angles, and the newly shocked gas at increasing angles can become 
dominant in providing the observed flux. This second case introduces additional complexities and 
will not be discussed here, since even the simpler case first mentioned above will serve to illustrate 
the point that a great variety can be expected in the temporal behavior of afterglows. 

For the conditions where the afterglow is dominated by the newly shocked gas at increasing 
angles, the observer-frame peak frequency of the synchrotron radiation spectrum from the blast 
waves (11) coming from increasingly larger 6 at increasing times t is 

'^m OC Csfe n'^^T^ OC t-[12fe-43fc+./2)]/[8fc-,-2d/c] (^3) 

The observer-frame intensity at this peak frequency is I,^^ = F^/^, ~ F^(£'esj;,m/4'/rr^tfm) oc 
S,]^^£,e'^esy,m t-'^('=-J'/2)/{8*:-i-2'ifc)^ where as in §^ the synchrotron efficiency is 
esy,m ~ {^'ex/'t'sy^mY OC (C^^e'^^~'^F^)" if synchrotron and adiabatic cooling are the two 
most important energy loss mechanisms (or its generalization if IC or other effects need to be 
included). We have then esy^m oc i-'i{4'c+j(i-c())/(8A:-j-2c(fc)^ where a = 1(0) if the peak electrons in 
the deceleration blast wave at the angle corresponding to detector time t are adiabatic (radiative) . 
(In this model the dominant radiation is produced at the initial deceleration blast wave for that 9, 
so a does not enter in the dynamics, only in the radiative efficiency of the initial blast wave). The 
flux observed at Vm from the deceleration blasts at increasing 6 is F^^ ~ t^V'^I^^, or 

i'um OC Vm 



OC 



^[Ak~2j^d(k-j/2)^a{Ak+j{l-d})\/[?,k-j-2dk\ ^ (^^^^ 



which scales with Ce'^^^e"^- F)epending on the normalization of eqs (p!o|), the flux ( [l^ ) from 
increasing 9 values can dominate the flux given by equations (Q) or (^). (In other cases, one can 
approximate the evolution as the superposition of isotropic blast waves from individual 6, which 
could in some cases be dominated by that of the central jet region). At a fixed detector frequency 
vd, the observed flux corresponding to eq. (|l3) is 



Fd ^ Fy^iuD/T^mT OC t[fc(4+12a)-2j-4A:(l+3a)-(i/2){l-Q)}-a{4fc+j(l-d)}]/[8fe-j-2dfe] ^ ^^^-^ 

which scales with ^^"~^~"^g~^~^". For characteristic spectra with a positive (negative) below 
(above) the break, this leads to detected fluxes which initially rise in time, and then decay. 
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However, a variety of behaviors are possible, including some where the flux after the break passes 
through the detector window continues to grow at a slower rate, or saturates. Note that the 
scaling of eq. (|14|) with allows both for Fy^ to decrease or to increase as Vm decreases in time, 
both in the radiative and adiabatic cases, depending on the values of j and k which characterize 
the angular dependence of E and F. 



5. Discussion 

5.1 Dynamics, Cooling and Decay Law.- In the simplest model where the GRB and the afterglow 
both arise from an external shock (case al of [Meszaros &: Rees, 1997a ), to zeroth order the 



GRB 7-ray flux should lie near the backward extrapolation of the afterglow, provided the basic 
conditions have not changed and the same radiation mechanisms are responsible for both. This is 
clearly a rough approximation, since it is likely that the GRB is initially radiatively efficient, and 
becomes radiatively inefficient at some later stage. The observations of GRB 970228, especially 
following the HST observations of September 1997 (Pruchter, et al, 1997) indicate an optical flux 
decaying as an approximately constant power law in time. Superposed on this overall long-term 
behavior, there may be shorter timescale variations (e.g. Galama, et al, 1997), which could be 
either a difference in calibration, or a real wiggle in the decay. 

At late times (e.g. six months in GRB 970228) the remnant is most likely adiabatic. One 
question that can arise is whether a simple external shock afterglow model whose dynamics is 
manifestly "adiabatic" can be radiatively efficient enough to produce the initial relatively high 
X-ray and optical afterglow luminosity. For an afterglow luminosity Lx,o ^ L-y, as observed, from 
our discussion in ^ this is not a problem. At each radius, the electrons can radiate up to half of 
the total newly shocked proton energy randomized in the shock transition, as long as the electron 
cooling time is shorter than the expansion time. The electrons can be radiatively efficient even 
when the dynamics of the remnant (i.e. the shell of hot protons and magnetic fields behind the 
shock, which provide most of the mass and inertia) follows an "adiabatic" law F oc r~^/^ for a 
homogeneous medium. Arguments were presented ^ why the latter behavior may be more likely 
than a faster evolution with F oc r~^, a conclusion also supported by comparison with observations 
relating to the size of the GRB 970228 remnant ( |Waxman, et al, 1997cD . Note also that, from 
equation (0), for expansion in medium whose density decreases, F could drop even more shallowly 



than the above for a homogeneous medium. As shown by P?avani, 1997|jWaxman, 1997a| ; [Wijers 



Rees &: Meszaros , 1997 ; Reichart, 1997| and others, in GRB 970228 the initial 7-ray flux and the 



overall X/0 afterglow behavior are in good agreement with a simple external shock where the 
observations started after Vm < t^d-, without any substantial changes of slope during the observed 
decay phase. However, as discussed below, there are various mechanisms capable of producing 
changes in the decay slope, which could be responsible for some of the reported departures from a 
simple power law behavior. 



5.2 Intensity Offsets Between Afterglows and Main Burst .- In the case of GRB 970228, a 
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backward extrapolation of the afterglow flux shows some hints of undershooting the 7-ray flux. If 
this were real, a slight undershooting could be due to the GRB being radiative initially (this is 
expected especially during the initial deceleration shock), with the afterglow dynamics becoming 
adiabatic soon afterwards. This would lead to a flattening of the spectral slope (e.g. second line, 
last two columns of Table 2). 



Another possibility ( [Meszaros &: Rees, 1997a ; Wijers, Rees & Meszaros , 1997; Katz & Piran, 



1997) is that the 7-rays could have originated in an internal shock. Internal shocks leave essentially 
no afterglow, yet they should be followed (eventually) by external shocks. An internal shock leaves 
unused anywhere from ~ 20% to ^ 90% of the total kinetic energy ( |Rees fc: Meszaros , 1994 ; 



Kobayashi, Sari &: Piran, 1997 ). The leftover energy is liberated in the external shock, most of 
whose initial radiation can come out at GeV energies because initially inverse Compton (IC) losses 
dominate over synchrotron losses; the initial synchrotron MeV radiation would then be typically 



low below the BATSE threshold (^ eszaros Sz Rees, 1994[ ), but later on it becomes dominant over 



IC ( Waxman, 1997b| ). Thus at later stages the afterglow from the synchrotron peak can have 



a flux level whose back extrapolation might overshoot the 7-ray flux. Since external shocks are 
generally smoother (at most 3-5 pulses, panaitescu &: Meszaros , 1997 ), while internal shocks may 



be very variable ( Rees &: Meszaros , 199^ ; |Kobayashi, Sari fc Piran, 1997| ), a natural conclusion 



(also reached independently by Piran, 1997) is that since afterglows appear to arise from external 
shocks, a burst where the gamma-ray light curve (or at least the last gamma-ray pulse of the 
light curve) is relatively smooth has a better chance of leaving behind a visible afterglow at lower 
frequencies. 

A reason for offsets may also be (§^ that the relativistic outflow has an angle dependence 
in either the energy, the bulk Lorentz factor, or both. The GRB itself may be due, for instance, 
to a high r ejecta which shocks first, while the afterglow could be dominated by slower material 
ejected at larger angles relative to the observer, which shock later and produce softer radiation, 
but which could carry a substantial or even larger fraction of the total energy. If detected this 
would generally be as an upward offset of the afterglow relative to the GRB, since otherwise the 
afterglow would be dominated by the evolution of the same material which gave rise to the GRB. 

5.3 Rise and Decay, Late Rises, Bumps .- An initial afterglow flux rise followed by a decay is a 
direct consequence of the simplest afterglow model (al), and is true also of any generic peaked 
spectrum from an expanding cloud where the peak energy decreases in time faster than the peak 



flux. Estimates for expanding clouds were made by Paczyhski Sz Rhoads, 1993 and Katz, 1994b 



based on simplified radiation models. In the more detailed model (al) of Meszaros &: Rees 



1997a , if the slope a at frequencies below Um is positive, oc u"', the initial rise is 



for ij£) < in a homogeneous medium ( [Wijers, Rees fc Meszaros , 1997 ). For an "average" 



GRB spectrum, a ~ below the break ( [Band, et.al., 1993 ), which implies Fy^ initially constant 



( Meszaros fc Rees, 1997aD . However, a is only the average value; there are many GRB with 



a > below the break (this is also the case for an ideal synchrotron spectrum a ~ 1/3 below the 



break, e.g. Meszaros , Rees fc Papathanassiou, 1994 ), and in such cases one obtains an initial 
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power law increase in Fq. After the frequency Vm of the peak has dropped below the observing 
frequency vd, if the spectrum above the peak is oc with /? < 0, the flux observed at vd 
starts to decay cx This occurs after a time tD/t-^ ~ (v^jvDf'^'^ (foi' adiabatic remnant 

r oc which for observations in the R-band is topt ~ lO^t-y. A maximum at 1.5 days such as 

in GRB 970508 is therefore compatible with the observed ~ 5 x 10^ s. 

However, there are other mechanisms which can give rise to optical fluxes rising at times 
which could be even later than the above. One is an anisotropic outflow (§§), where the GRB and 
X-rays come from material close to the axis oriented near the l.o.s., and further off-axis material 
with a slower F starts to decelerate after 1.5 days, say, or it has been slowed down enough that 
its light cone includes the l.o.s. after 1.5 days. Since beaming does not change the slope, if the 
spectrum is a pure power law then its slope would remain constant. While a constant spectral 
slope is the simplest assumption, it is not a necessary one. It is conceivable, for instance, that the 
slope of accelerated electrons (and consequently of the synchrotron spectrum) could change as the 
bulk Lorentz factor and the shock strength changes, or it might depend on other effects associated 
with the angular anisotropy of the outflow. An interesting result is the indication that the optical 
light curve of GRB 970508 may have been steady or even decreasing ( Pedersen, 1997 ) before the 



1.5 day rise phase preceding the maximum. This may be simply explained in an anisotropic model 
such as described in §^ where the indices j and k change to give a light curve transition in this 
sense, or even more simply, by a bimodal model where one has a central jet associated with the 
7-ray event, whose tail is just seen to decay, followed by the emission of a slower F outflow over 
much wider angles outside the jet which is responsible for the main part of the afterglow. Another 
possibility for a late rise in the optical would be if a lower F shell catches up with the main shock 
front with a comparable energy after t ~ 1.5 days, so we see then the emission from the onset of 
deceleration of this late shell. Very prolonged optical decays with a shallow power law are possible 
in anisotropic models such as in §^, e.g. Table 2, lines 6 or 7. 

An alternative explanation for a late turn-on of an afterglow may be that the GRB occurs 
inside a very low density cavity inflated by the pulsar activity of one of the neutron stars in 
the progenitor binary. The shock, at least over a range of directions, would not arise until the 
ejecta hits the wall of the cavity, and this could take a time of order weeks, the ensuing shock 
being spread over a dynamic time scale sufficiently short to produce a large flux per unit time. 
A characteristic feature of pulsar cavities is that they are usually asymmetric and irregular in 
shape, often being elongated due to the proper motion of the energizing source. One would 
naturally expect a wide variety of time histories for the afterglows arising from the impact of a 
(possibly anisotropic) ejecta upon an irregularly shaped cavity wall whose dimension (depending 
on direction respect to the line of sight) may vary considerably. 

5.^ Peak Flux Level Evolution.- In some observed afterglows the flux level F^ at lower energies 
is, at least initially, significant relative to that of the maximum gamma-ray flux (even if lyFi, is 
smaller). One prediction of the simplest model is that the maximum value of the afterglow flux 
in every band, F^^^, is a constant. An interesting case is that of GRB 970508, where the ratio 
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of maximum Fx to maximum is ~ 1/2, while the ratio of Fq to F^ is ~ 1/10. Considering 
the drastic simphfications involved in the simplest model (homogeneous medium, constant 
equipartition fields, etc), this order of magnitude agreement is perhaps encouraging. However, 
-Fi/^ ~ constant is clearly an approximation which need not always hold for less simple models. 
In Table 1, the top four lines of the last two columns show that for isotropic but inhomogeneous 
outflows one can expect F,^^ to either increase or decrease with as the later decreases, while 
the previous columns of the top four lines in the same Table show that either of these behaviors 
may also arise as a result of an anisotropic outflow. Thus, the ratio of -Fjy,,^ at 7-rays and lower 
wavelengths could be either larger or smaller than one, simply on this basis. In a simple bimodal 
anisotropic model, one can simply have more or less energy in the large angle slower outflow seen 
later than in the early and harder axial outflow. 

A related question is the observability of radio fluxes of order mJy around 10^'^ Hz, as reported 



for GRB 970508 (Frail, et al, 1997). Radio fluxes of this magnitude arise naturally in a simple 
isotropic homogeneous fireball model after a week or so, since the self-absorption frequencies 
(overestimated by 10^ in Meszaros &: Rees, 1997a| ) are in the range of 10^^ Hz initially and drop 



to lO^'^ Hz in about a week, and the flux level is well below the brightness temperature limits 
for incoherent synchrotron radiation. What is more interesting is the relatively large value of the 
radio flux (~ mJy) relative to the O flux ('^ 50/iJy). In this afterglow therefore -F^m first decreases 
with between 7 and O energies, but then increases with decreasing between O and the 
R(adio) band. The simplest explanation may be in terms of a jet-wind two component model. 
The decrease between 7 and O could be due to expansion of a jet into an inhomogeneous medium 
(e.g. Table 1, fourth line), and the increase between O and R could be due to a surrounding low 



r wind at larger angles, which shocks at later times, as suggested in Wijers, Rees & Meszaros 



1997. The slow wind need not have an angular dependence; a growth and decay of the flux is 
approximated also by behaviors represented in Table 2 in the last two columns for expansion into 
either a homogeneous or an inhomogeneous external medium. A wind with F ~ 3 — 10 can match 
the delayed emergence (~ week) of the radio from the wind blast wave. The tables give only 
illustrative values for selected spectral and density exponents, which can be easily changed to fit 
a particular observed rate of growth and decay. Actually, the radio flux of GRB 970508 at lO^'^ 
Hz (Frail et al 1997) at first increased and then appears to flatten, except for decaying oscillations 
which could be due to scintillation ( poodman, J., 19*97 , Waxman, et al, 1997(j ), followed by a slow 
decline. 

5.5 Burstless Afterglows and Afterglowless Bursts.- An interesting consequence of anisotropic 
models (Meszaros & Rees, 1997b, Rhoads, 1997| ) is that there could be a large fraction of detectable 



afterglows for which no 7-ray event is detected. If the observer lies off-axis to the jet, then the 
detected "afterglow" can be approximated by the isotropic model calculated for an E{9obs), ^{(^obs) 
corresponding to the offset angle 9obs of the observer to the jet axis. As F drops after the 
deceleration shock, the causal angle includes an increasing amount of the solid angle towards the 
jet as well as towards the equator, and depending on the values of j and k the observed flux would 
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generally decrease in time. However for some choices of j, k an increase in the light curve might 
be possible, depending on the normalization. It could be that gravitational energy is converted 
more efficiently into kinetic energy of expansion at large angles, where the opacity is larger (c.f. 
Paczyhski, 1997). After T{9obs) has dropped to the point where the central jet portion ^ ^ is 



detectable, the late stages of the jet emission would become visible, at a later stage when it is 
bright only at wavelengths longer than 7-rays. If the jet contained substantially more energy than 
the off-axis regions so that it dominates the flux even after expanding for a longer time than the 
initially observed off-axis region, one would expect an additional increase or flattening of the light 
curve at this point. Details would be further complicated by contributions from the equator and 
the back side of an opposite jet, if 9obs '/r/4. The statistics of afterglows not detected in 7-rays 
can be calculated from equations (|ll|, |l4|). 



The converse question is why some bursts (e.g. GRB 970111) have been detected in 7-rays 
but not in X or O, even though it was in the field of view of Beppo-SAX, which would have been 
expected to detect it if the X to 7-ray ratio had been comparable to GRB 970228 (a weak X-ray 
afterglow may have been detected, Costa 1997bD . One reason may be if the 7-ray emission is 



due to internal shocks (which leave essentially no afterglows, Meszaros &: Rees, 1997a| ), and the 



environment has a very low density, in which case the external shock can occur at much larger 
radii and over a much longer time scale than in usual afterglows and the X-ray intensity is below 
threshold for triggering. This may be the case for GRB arising from compact binaries which are 
ejected from the host galaxy into an an external environment which is much less dense that the 
ISM assumed for usual models. Another possibility for an unusually low density environment, 
made up only of very high energy but extremely low density electrons, is if the GRB goes off inside 
a pulsar cavity inflated by one of the neutron stars in the precursor binary. Such cavities can be 
as large as fractions of a parsec or more, giving rise to a deceleration shock months after the GRB 
with a consequently much lower brightness that could avoid triggering and detection. 

The lack of an afterglow in some bursts may also be due to occurrence in an unusually high 
density environment (e,g, a star-forming region, or the inner kiloparsecs of a late type spiral, 
where failed supernova or hypernova progenitors may reside, e.g. Paczyhski, 1997] ). This could 



lead to a more rapid onset of the deceleration leading to the X-ray phase, and it would also imply 
an increased neutral gas column density and optical depth in front of the source. A special case 
is that of GRB 970828, where X rays have been observed, but no optical radiation down to faint 
levels (Groot et al. 1997). The presence of a significant column density of absorbing material 
has been inferred from the low energy turnover of the X-ray spectrum ( Murakami et al., 1997] ), 



and the corresponding dust absorption may in fact be sufficient to cause the absence of optical 
emission (Wijers & Paczyhski, private communications). The difference between the low density 
and high density environments cases could be tested if future observations of afterglows reveal a 
correlation with the degree of galaxy clustering or with individual galaxies. 

In conclusion, the absence of detected afterglows in many bursts is not surprising, while there 
may be detected afterglows also in some cases where a corresponding gamma-ray burst has not 
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been detected; and when afterglows are detected, a wide diversity of behaviors may be the rule, 
rather than the exception. 

This research has been supported by NASA NAG5-2857, the Royal Society and NATO 
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Table 1: Exponents of the power law dependence of the synchrotron peak flux Fi,^ as a function 
of the time- varying synchrotron peak (top), and exponents of the time dependence of on 
observer (detector) time t (bottom). The first column indicates whether the electrons are in the 
radiative (a = 0) or adiabatic (a = 1) regime, and the second column indicates the value of the 
exponent of the external medium density dependence on radius n oc r~'^, d = being homogeneous. 
Columns 4 through 11 give the exponents for the anisotropic model E oc , T oc 9^^ of §^ and 
various values of j and k. The last two columns on the right gives the corresponding exponents for 
the isotropic models of the first being for the strong electron-proton coupling a = A case, and 
the second for the weak coupling case A = 1. For A = 0(1) the remnant as a whole is dynamically 
radiative (adiabatic). 
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Table 2: Exponents of the power law dependence of the synchrotron flux in a given detector band 
as a function of observer (detector) time t. The first column indicates whether the electrons 
are radiative (a = 0) or adiabatic (a = 1), the second column gives the exponent of the external 
medium density dependence on radius n oc r~'^, d = being homogeneous, column 3 gives the 
value of the spectral index a of the spectrum, F^, oc z^", which may be 1/3 below the peak j/^, 
and -1 above (although these indices can vary around these representative values). As the peak 
passes down through the observation band i^d the index a changes from the positive value left of 
the peak to the negative to the right of the peak. Columns 5 through 12 give the time exponents 
of the Ff) cx in the detector band for the anisotropic model E oc 6^^ , T oc 0^*^ of and various 
values of j and k. The last two columns on the right give the corresponding exponents for the 
isotropic models of §^, the first being for the strong coupling case a = A and the second for the 
weak coupling A = 1. For A = 0(1) the remnant as a whole is dynamically radiative (adiabatic). 



